{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "a045577b",
   "metadata": {},
   "source": [
    "# 栅栏效应"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "b60eb4ce",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEGCAYAAAB1iW6ZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAABVjklEQVR4nO29eZhk2VXY+buxZ0RG7kvtldXd1Vq61epulVobFrIZUEsICSNsw+DBxmY0wmjYDcLM8IlvjD2AhwEPiyxkMAIhsUhgISS0YGHJaqTuavVevVV37VW5RWRmbBn7nT/eexEvIt9yX2bki4rs+/u+qoyMiPfi5I1zzz3n3HPvFVJKNBqNRjP6RIYtgEaj0WgGgzboGo1Gc0DQBl2j0WgOCNqgazQazQFBG3SNRqM5IMSG9cFzc3NyaWlpWB+v0Wg0I8nDDz+8LqWcd3ptaAZ9aWmJs2fPDuvjNRqNZiQRQlxye02nXDQajeaAoA26RqPRHBBG3qCvFWt8928/wJ8/cjXwtfVmm/d/4nF+7s+f2AfJ1PjsEzf4J//p77iUKwe+9vxqiX/0wQf4b8+s7INkavzUnz7G//kXT9JqB19x/PEHL/OPP/h3bFUa+yCZP9c2t/mHv/VVPvvEjcDXVhstfuKPH+UX/+rcPkimxicevsr3ffhrXN/cDnztuesFvuu3vsoDL6zvg2T+SCn54T/6Bv/20+fYzWr13/vqBb73Q1+jUm/ug3T+XFgv852/+VW+9OzqUD7fjZE36H/28FXOXtrgx//4MaqNVqBrP39umY8/dIWPfv0yl3OVfZLQHSklv/K5Z/n6hTy/+oXnAl//S3/9DA9d3ODXvvj8Pkjnz7PLRf7s4av8wdcu8bcBFbtYbfD+Tz7Bgxfz/MWj1/ZJQm8+9vXLPHJ5kx/66DcCD0j/9dFrfPKRa/zOVy6wWqzuk4TutNuSf//ZZ/jq+Ry/8aXzga//t391jm9c3uQ3d3HtIHj0yiZ/9fgNPvw/LvD1C/lA164Wq/zCX57j717M8ZknlvdJQm9+/4GLPHplkx/52CO7GpD2i5E36J97qvuFPnplM9C1X3sx13n8+XPhK8aV/DYvrpc7sgRRjHZb8qDZER6/ujUUo/J5W9vb21KFsxc3Oo/t32GY2D/36RuFQNd+7cWuEfpvT4fvpT2zXGS9VDNlCdb2tWaLhy8Z7f/ACzlKtfC93M891Y0qg8r/dVvb//WTw9EdS/eL1SYXAzqD+xlVjLRBb7clz68U+a57jiIEPBRwpH/owgbffPs8c+MJnlsp7pOU7jxrfuZ3v+YYK4UaV/LqofNzq0W2thu8+95jADy/UtoXGb14dqXI8Zkx7lua4UGbgVbhwYt54lHBd7z6yFDavtZscWG93Gm/BwPqzoMX8rztzkOMxaOd7zFMnrPpzotr5Y5xV+HJa1vUmm3efe8xpDRSd2Hz7HKBlx/KcseRCR66GLDfXsyTSUR56x2LQ9GdQrXB9a0q3/0aQ3eC2B0pJff94t/w7z/z9L7INtIG/drmNuV6i9csTXN0aozza+qK2Wy1eXG9xCsOT3B6IctzQzCIljK+89VHAHh+VV05LQP+nfcc6blXmDy/UuL2hSyvOJzlhdVSoAjj/GqJU3MZXn1skvVSnVwAgzQILqyXabYlb759jul0nOcDGLVKvcm1zW3uODLB6cXxoQymz60UiUUEb3/VISCYUd6hO8vDGJBKnF7M8orDE4EHlPOrJW4/ZFx7ZaPCdj1YqnWvPG/2tW975SLJWCRQv10v1SnVmhyaTO2LbCNt0C1FuH0xy9JsJlDoc2OrSqMlOTWX5vbFcc4HNEiD4PxqiSOTKe48OgkQSH5rEvU1J6eZSsdDH5CsAfH0YpaluQylWpNcua58/aVcmZOzGU4vZgFCl98yaqcXspyczQSalL5kfk9LcxnTGRjCYGoOiKcXsqZM6vJfzFWIRwWvOzVLMhYJXf5yzRgQb18YZ2k2zUqhFigNcSlXYWk2w+2L2aFEGJbuvPzQBCdn07vqt0tzmX2RbaQN+uW80ZAnZ9MszaW5uK6u1BfM956czXQM0kbI1RaX8xVOzmaYTsfJpmIB5a+wOJEknYixNJvh6ka4k7rLBWNAXJpNszRrKKeq/O225FKuwqm5DKfMa6+ELL+lO0tzaU7NZQK1vfXepdkMp+bSrBZrgSfk98oVU3cOT6aIRwUX1tXb7+J6meMzaRKxCCdm0qG3/dUNI7V4ci7TMWyXFI1itdHi+tY2S7OZjt4NQ3diEcHR6TFOzgbTnQs23dkPBmbQhRDHhRBfEkI8LYR4Sgjxo4O6txs3tqrEo4K5TJKl2Qxb2w02FL1Ea6Q8NWd0CuN+wcu/9sLyVpXDUymEEIZRCeQlljtKcXgyxY2tcCdFl83POzw11umUFxQVe7lQpdZsc3I2zcJEsud+YbG8VWVyLE46EePkbJrrW1Vlo2x5ZCdn0xyaHANgpRCu/De2qhyZShGLRjg+E8yZuZgrdwbSQ5Op0Nve6mdHJlOBnYEr+QpSGgNxt9+GrzuLEymiEaPfXspXaCtWSV3KVYhGBMemx/ZFtkF66E3gJ6WUrwBeD/ywEOKVA7z/Dpa3tlmcSBGxNdA1xZrcq5vbJKIRFrLJTqcMU7FbbclKodpRymPTY8qyg/F3HptOA8PqlKZBn0xxdMpov+ubajJYf+ex6TSpeJTZTCL0Tnljy9726c5zKlzbrDCVjpNNxYdiVCr1JlvbjU4e9th0musBnBFDd4zvbJjOwKHJ1K76LRj9ZSodJxmLsByyI3Zjq2pr+zHqzTbrZbU5oGub2xyaSBGP7k9yZGB3lVLekFJ+w3xcBJ4Gjg7q/k7c2KpyxDTGixNGA6t6SitbVRYmkgghOp3yeoiKvV6q0WzLzmCyOJFiRfHzW23JarHGoundHpkco1RrUqiGlzKyd8pELMJsJsGyYtt3rjW/s8NTqdA75XJhu9MpLTlUB8XlrVrnmkNDiO7sgynAoYmksuyVepNitcmC1faTY6yVajRa7f0R1oEbW1WEMHR+cswwykH6LRjXCiE4MjUW/oBU6Br0haxpd7bUDPqybTDYD/ZlmBBCLAH3AF/ve/49QoizQoiza2tre/4ce8NaP5WNSqHa6ZRz40miERGqUbGU8IjNqJTrLYoKRjlXqtFqy51/e4iKfWOrSiYRJZs09ndbnEipd8pCr0E/NDGETrlV5bA5mB6aTPbI5cdKodpxIKy/IUz5uwPiWEeG9VKNpoJR3jGYTqaQElaL4VUZLW9VmR9PEo9GEEIYEWZB0SCa35FlSA9NhBudSim5sbXd7bcB7c6Kze7sBwM36EKIceATwI9JKXtWa0gpPySlPCOlPDM/77j7YyCMjmV0xvnxJBGBspe7UqixaH4Z0YhgIZtkWXGUHQTLNk8DuoqhYlQs5em/NkzFtoyaEKIjg7qHWyUVjzAxFjOvTYaag2602qyX6h3dsdpxN85AJhkjm4op690g6OqOKf9kiraENYXST+tvtHRmcQi6s2wbECFYdLpSqDI3niARM0yXMRiEJ3thu0m10d4xoKvIIKXc8bcPmoEadCFEHMOYf1RK+clB3ruf7XqLaqPNTMZQ6lg0wtx4Ur1ht3pHyul0gs2KetndXsmZObe58T6jojCo9HtZ0+kEABshyr9eqnVkh2AeumUQrcFgJp1ga7uxq/1gdkPenDi35M+m4mQSUSWjZgwGXWcAYCaTCLVCqqM7WUP+ICmjlT5nYMbUnbB1f2480fn90IS6UbYmJC2Mfhte26/39du58YSyI1msNanUW52IcD8YZJWLAP4z8LSU8lcHdV838qYCzmTinecOTaZYUQjdirUm241Wx8Mx7pPo3DMM8iVLfkOxrU6pYhRXzPC40ynNe6hW+AyCfLnObF+nzJXr1Jv+Yf9qodbJ4QJMZxK0JRS2w+mYuZJl0LvyL06mlLZPWC/VkJIe3ZlOJ0IdTHPlOolopCfdBSjpvvUeS35Ld/Jh6k6pzqzNGTD6bVVpHchKodZj0GcycUq1JrVmOGWjVjtZuh+LRpjPqkWYq32D6X4wSA/9TcD/AvwDIcSj5r+3D/D+PVjGy/JOAWYziY734kW+1OuhgWFUwjSIuXKdbCrWCR0tBQkiv3XN5FgcISAfqpdY7xgDuywqhi1XrjE/3juYAqENqFYbW9EdwFwmyXpJQXYH3ZnJJEI1iLmSMZhaEY4li5LulOskYxHGzcFgOhNudCelZL1cZzbT229rzTZlhRWf+XK9ZyC25A/LS7dWNPfofiaptKjOSXcGzSCrXP6HlFJIKe+SUt5t/vvMoO7fj9WBeo1KsmPsvLAaf9p27Uw6Hmqn7E9ZjCcN466iGPlyjWwq1il9ikYEU2Px0AakVluyUen1sqxOllNo/3y5zrQtsuqkjEKSP9c3IFqPVb5/J72bTofsDJRqPbJ3BkTFAWkm0x0MMokoiWiEfDkcg1iqNak3231tb+iRn/xSSlN37P023Ahj3cEoz44nFPvtTid00IzsSlHLo5jq89DXy3Xf0K0TNmV6R/pCtalUKTAI8n1eihDCiDAUByT7tRBu2L9RqSNlb8rC8nb9vMRWW7K53ejxjqdD7pRW55vL9HrZKvvJOBn0mUw81Bx6vlzvab9ELEI2FVN2Bmb69G4qHQ8th97tezaDaMrjV8tdqjWpt9o7+i2E5ww4GeVZRd2xvh/7YDZoRtagO3voCeoKoduGY6c0Q7cQ87gzfUZZNXTfqOy8djoTnkHPlZzbHvyN8qY5GPR2SsNbD0/+GrGI6FTZgOElbm43fAd0N2dgu9EKbZOo9VKdub7v30g3Khj0SmPXejcILA93ZtxBd3ycmQ0zipjJDDFdV6oxYUuVgpkZUOm32kN3Z6PSQAgjf2zR8RJ9Rsucg5cwFXbYX95plGfHk2ojfanXQwNDScIKmzs56D4vBfDNQzt7uFYeNxz58+U6U+lu2gEM+aX0lyFfrhONCCZSDimjEL3c6V3qTr5cG2p01/n+0zu/f7/oznrdLv9U2nIGwtJ9Z0esUvcf0HPlOtlk72AwaEbXoJfrTI3FiUZsnbIzsehnVGqk4hHGEtHOc2Hn4grVBpPpeM9zyl5Wud5T3QNm2B+W7NvGznh2+SdScWIRQd63U+406GPxKMlYJDz5q42OIbBQjTBy5TrT6TgRm96FmTKqN9tsN1pMjfV//4pzAKWdg0GYHrpVyWRv/9lOuk7NGbDLH/b8S6HaZLLPw55TLGjIl+s9kcl+MLIGPV+p7whdrJHbLw9t5KD7PNwQw/5qo0W92e6JLgClHLqU0ky5OHjoFf/5g0FgdUq7/JGIYFpBficPXQhhRhjhdMqt7caOtu94iT5ebn8O2n5tGLpjbe/Q7wzMjSd8o6Nqo0W53trhoU+lw5sD2HLQnbFElHQiqtRvoddDj0eN+YPh6o6VGfDX/X7dGTQja9A3XMJOwNdL3HAJmyCc0M0yiPawHYy8ol8utlhr0mjJnWFzJtHx3vYby6hMOA1Iil6Wk/xhhf2F7SYTqVjPc93SP/887k7dCS/sd9Uds/28dv3b6Kzd6HUGZjLGojrVHQP3gqU7VtmkhUqVkdPcl/V7WLpT3G7s0B3V6C5frvekmvaDkTXo+bK7h66Sx+0fDMIMmy0vxckggnfoZk0c7QibQ5ZfCBhPBO+U1utT6Z1GMUwvq7/tVRfY5Bw89DDD/q7u9Lb9TCZJqy09N2jrTmb3/u3TaXNhVwibu21tNxhPxoj17TY4o1DLnS/XScQipG2pUiD06M6932oPfdcYaYfehk3Fo2QUQ7d+DzEVjzIWj4bSKTth8w7F8A/dnMJOsJdvheMlTqR688iGTP4Tc/m+BVUWxsRcOGF/obozbJ5OJxBCJeWys1Na9wrDqBSqzZ7PtLDyuF7OTDfdtdNDt7++nxS2mztkB5hTKP2z+q19MhtgOh0PxUOXUjrqjpUZ8JLfqqHXOXQHjDxyw7H8Z2Y84ZtycRspw1r+3/Gy+kK3GYXQzSkHDYZSQzjlW4aXEtvx/IxCysVpMIXwKi3abdkZkOxEI0Ye30t+pxp6MJZ/T46FY1S2PFIusDvd6VaKhCN/NuWsO0opCyfdySRCcWS2Gy0aLbmj7TOJKIlYxFN+pxr6/WAkDXqlbkwq9qcdwD90qzZaVOotF8WIh7KEuFMl0u9lmYbC6wR3tzxidwl0OF6ik5c1m0lQrHrvq+E0fwGG/FsKdeB7pVxv0pY72x6sxUXu7edUQ2+/Nswcunt055Guc4nuZsKM7hw8XDCcmVzJe1Lf1RELyRlw67fWokCv6Mhq2/2sQYcRNejd1VpuoVtwL8W4n1rZ4F5xm1RU8dCdyv6gm0NXWWm6V7YcPFzohp5ehsGpjheMrRek3P+FXVbKwinCmPXxEp3K5iym03HfyHAQuE5IK5Ts5st1ImKnQQpz/qjgkIMGw5mpt9qUau6HRXt56JV6a9/Pde22vYPu+GQGOjX0OuWyk27pk7On5Dmp6GHQJ8fiFAMalE88fJWf/eTjgSoEtirOYXMmYdRje3fKGkmHiSGrk2wFkL/ZavOTf/IYn378uvI1YHRKNw8XvCMMp7I/6JbhBdlx8bmVIj/4+2cDHd1ntb1jhDGe8Fx+7jZ/Yd3P8uBU+cOvXeIDn3oqUKnp1naDRCxCKr5zYhD851+m04kdcx9W2wfRnVqzxY987BG+eG5F+Rrw1x0/Z8yt31r3VuWJq1u85yNnlXbYtHAqubTwywy4zV8MmpE06F4j5YxZaeHWSfw6ZRClXi1U+ck/fYyPPXiFTz5yTfm6QrXBWDy6Y2JQCOEb9rtNDEUjgmwyFqhS4WMPXuYT37jK+/7okUCTwe4euren15kYclDqyV0MSP/6Tx/ji0+v8H/95Tnlazq64yC/Xx7XzxkIIvvlXIX/4y+e5L88cJHPPbWsfJ1RcrlTdms/Fy8v0W0wHU/EiIhgVS6/99WLfOqx6/zgR85S9vCq+3HTnRmfCKPWbFGqNV37rXVvVX7k44/w+XMr/PJfP6t8jVvJKPivIelE1jrlshPLE3Jr2EZLuoZubjloMLzcQrWh7DH93Yu5zuO/fXZV6Rpwn1S05PLKB254zJRPBPQS//bZ7jGAZy9tKF/ntMoV/CfmSi419ND9Lq2UiB9b2w0eu7oFwFeeX1M+HMOtZBQM78krj++Wg7buF8Qg/o/z653H9u/Bj4KH7sxmEp5bKDvV0IOxKCybigfycL/0TFffH7uyqXRNs2Xss+Q2/wLuuuO0j4uF9V2qtv/yVpUL62XAaHvV/u6tO/79FtBVLk64lf1B9wt3Uwy3HLR1v0ZLKi/OeeTyJmPxKPffcYhHLm8qXQPupVuWXH6hm9vEykQAL1FKySNXNvmOVx8hFhE8clnNoNeaxklR/RU64N8pvXLQQb2sx69uAvBd9x6lXG/x/GpR6Tq3SUXo7ufilsd3q6G37lfYbiin3h65vMFMJsE33z4fTHdcJhXBijC887huddBBIoxmq83jV7f4h/cYZ8A/omjQix7zF36HtHT3sN/5twfVnUevGLr+XfccZb1U4+qGWsrOS3es/Vzc8vhWDX2mL1U6aEbToPuEPuAeuuXLtR2bK1l0vERFL/fxq5u86ugkZ5amuba5rbQ5EriHnWBNzHnncd1KnyZS6imXqxvb5Mt1XndqhlccnuBx09v1w22m3/h8Y28dv8HUzcM17q9q0A15/9kblozfr6jJ71b2B/4RhlsNvXW/tjSqaFR4/OoWrz42yT0npnhutUhF8Tov3ZnJJHeVgwbDyKpGR+fXSmw3Wrz59jlumcsoe+heOWi//Vy8ctCWc6Habx+7ukUsIvi+158AUNb9LfP+TmWXfnbHLVU6aEbToFebCOHcsJb357YVZ77ccJwYguAj/cVchVsXxrl1ftz8vawov7uXNZ1JeG4juuGSg4aul6iCJettC+PcOp/phKB+uFVZgLmfSzruqtRe6a6gbX9hvcxCNskdRyaIRQQXlNveXXf8Jub697C3E0T+dltyMVc2234cKeFyvqImv8ukInivtrVq6L3kV9Z7U1dum89yy/x4IL0H58F0LBElFY+4OjN+8xegrjuXcmVOzKR5+aEJIFi/zSSinYNl7PjZnQ2PyHqQjKZBN5cPOxnlTtjvks9y2j7UwgoFVbzcYrVBvlzn5Gya4zNpIECnrDqXboEhf9kldOtsruSZQ1dVakPWk7NpTsxmuLG1rXQeqFceEbzDfq90VzIWIRGNKEcYl3MVTs6miUUjHJseC2QQ3XTHb5MtPw/XuL+/l7harFFrtjkxm+HkbLrz9yjJX216zL8kzcNHdqZ9rBp6V/kD5NAt3Tkxm+bkbJrL+YpSHrrjoTvMv4D3UW5e8xfZVLDo7lKuwonZNJlkjLnxpHrbu5Rc2uVyszu5cn3fSxZhVA161Svs9A+bpx3ycND1HLYUFohYBuTETJrjM2MI0VV0P7YqXl6WWcvtoBh+R1gF8bKu5CskYhEWsylOzqRpS5TK/7zSXYb87qv2vLwsIUSgAelyvtIZSE/MZgJ1Sre2VwmbvQwiqHmJdt05OZPpec4LKaXjbn8WVkFA0aEgwGv+AoLpzuV8hcmxOJNjcU7Opqk22qwV/dONXsUMYOmOe791qqEHo8JnLB5Vkl9KyeVchROm7pycTXMpr+ahe7V91+64Rxj7vY8LjKpB3266jpRps5bbK4/bv3WuRaeeVcFLtAzIiZk0yViUwxMpJaPSbhsdzmlSEbzDfi+DCEZHKddbSqstL+UqHJ8eIxIRnDC9xEsKoWc3D+ouv9s6AOuA4v4a+o78YzElD7faaLFcqHaM4YmZMSXZwdsZ8E/XuU8qBqm0sGQ9OZNmMh1nIhVTcgYq9Rat9s6l5xZeZ4s6HerSL79ydJSvdCILa1C9pDAguW0sZjHtUTbqVkNvMako/2alQbHW7Bj0EzPpANGR19yX9z5MXsUMg2Q0DXp15xaWFn5nc7otPYdgE3OWN3tsegyAI1Nj3NjyX6RQrDWR0j1l4VXL3Qk7XVMuVsrI3yhe39rm6LSh1EenjL9hWUH+7kpLr5SLdw7abWJoIqXWKa12Pmq2/dGpNIVqU6ke2svLikcjTLjUckspzbI/H2dAQXeubxryHzHb/eh0Wkl3vCYVwbuW22v+AoyJxWqj7bltg8W1ze2Ozhwzf6rI71WdBt7bL3v1W1B3Bvr77dGpMVaKNaWy1y0PR3JiLOZaEOBVQz9oRtOge+SywH1v7e7mSu5KDd3ZbC9Wi8aKTUs557NJ1hSqXAo+OWiv48z8PPQgRmWlUOXQhGGcrAFCLWz2S7kYZ3M6dRC/3eZUw/6VgmE8Dk2kAKPtwXuFqoUR3Tk7A2Bt0LZTBr/NlYKs1F0pVpnNJDrVMsq64zEhDd1FK05pC78Diru6o6D7hRqLfW2vojtb2w1iEcFY3DlC83IGvNJdoK471spQu/yttlTaC8ZrDYB1SIvTfTo19DqH7kyx6rxazsKtlnvDZ2IoFjXqRFW8xJVClcWJVMfbnM8mlZUa3A2i16lLfqvNVPO4rbZkrdjtlMlYlKl0XHlASjosPbfLb5zN6Sy/V9ipmkO3DPqiOSAFNSp+uuPkofvloLPJGEKoRUerhSoLZtsDzI8nWVeR3WXLCLvsdlntdGvo3TxMNd0p1ZqUas2O7kyOxYlHhbIzMDkWd43QvGq5/XLQqtHdSsGQczcDktMunXbcMgNOZ6HuFyNp0L1WWoJ7LbefhwvqI/3yVrXjIYLRKbe2G74bBPmFnZNj7rXcVg2967VptTzueqlGW3aV2pJ/taBoEH2iI3D2Er0qjMDIywfx0BcnzU45HqBTepSMgnstt1cNPZirLZMxpQFp2RYdQdcZ8KsUcdsL3cJrg6582TigOBlzm79Q051OdDRpyC+EYH5c3Znx0h2/AWlQ/VaIriFXNegtc+7LW3ecIwy/YoZBMnIGvWnuyObtZSUdJ4a8Sp8sVL3E1WKNBVuntB77hf3dlIvzgORVy23U0O88WKIju+LCqK6HazPoAcJ+L6X2qhTxykGD5WU1fQ3bSqHGWDxK1jzGzGp7P/kbrTaVesvTqLjtuOiXg4YgEUZtR9vXW21fg+Q3qZhOxFxruf3SXROKpX8d3cnuRnfcc9DgbtBbbclmxX0NAATpt1VmM8lOLfmCokEv+qS7oLuPVD9+c1+DZOQMurVHi7dixB1ruTseuo9iq3qJ/Z0S/BXDa6WlxXTauXzLq8rCfk8/+bth504v0Q+n8zjtuHXKzsSQTw691ZaUPc5UNeSvsjiR7ITu0+kE0YjwjTC8lm5bWPMv/YOKVw29hYruNFtt1ku13pSLsu74yz+TTpB3KBtV8XDBX3esNu6Xf7WgMCnqcB6nHTdnYGu7QdsjVQrG/Fex1vTdesEYTLt6b50luzqAfjuTdj4gJ6ydFmEEDXq3ltXLqDjXcqvseGaUb3l7uJV6k0q91VEGgPlxQ8F9O6XKSO8RunnnoNUWRllRhF3+BcWw36tKBNw7pUrYqVpllCvVe2SPRozKJv+2d99LxMKtlls1XefX9nlzHmfeNrCpeonWvbNe0anLvty5svcBxaoVUpbuzNt1P5tSnJD2S3c513Jbv7vNX4ChO1LiWINvZ71U69GdTDJGJhFV77c+zsxmZefmbvlyHeFSQz9oBmbQhRD3CyGeFUKcF0K8f1D37UfVIMLOicX1Yg0h/BTDPw/qlLqxvCy/kX5ru0HE4YBlO7PjzrXc/Yasn7F4lHhU+HpZTsZpPptku9HyPGAAvFe5Qnfjqv6U13rRP+xUXQfg5G2qhP1+ZX/gXsu9XqyRirvX0INa6ZzTroFBdCebNMrj3OVPOjoDuVLNs+1VUy75ct3YC8k2KM5njRWefusf/HTHrd+umbrjpfuqzoDT9g2D0h2rfftPrlovGYOp1/c2KAZi0IUQUeA3gbcBrwS+VwjxykHcux+/sjlwr+VeK9WYSScc92KwUNkPxaniYXbcOGRYJWzOOhywbMfNQ18t1jqd3wkhhNIS7ny5TjoR7alUUQ37/Tx0a1/u/uhorWSE5Ase8quu1HUqYVNJGanojlst91qpxkI25bm5ksrEXK7jbXZlCJKu8zKI4FzL3WpLI81jy3v3k4obC/JUdMc4ULvbDvPZJFJ6n3jkt8oVupu77dQdM83joTuqKSNjpbiT7ninjPzKjcE93bhWrHr220EyKA/9PuC8lPJFKWUd+DjwrgHduwevwy0s3PbkWPMxiGAoVbHW9Fxo4JRPjUcjzKQTSiO9l+zGfXfWclfqRrmYfSLWCRWj4ujhKqSMpHQ+YLkfJ6Ni3der/bseuruX2zZrhvvlX1Aw6H770Fiyw84qHVXdUYkujM/p3suoPoko6Y7TpmJ2nJbP58t12pKB6I7Tbp/zCnnoaqPteMCynYh5UPdOg+ivO9399N3ldztPeCGbGojuuBv03jmT/WRQBv0ocMX2+1XzuR6EEO8RQpwVQpxdW1Pf1N/OiZkM/9s339JTMtiP2/maKp3SMipFD8WwwnHH0E0hj+uXS+ucr2kbkKyUxbxH2AmQVZgDcOyUlpfoYVRKNeOAZf8BaWceV6lTjlkLu9zbvlg1BlsnD329VPOcFPMrGQX38zXXijVPD9G6b6XeouGRenBKdwkhFHXH28O17ttfENBpex/dUVn+77RiU0V3/Cp0LJxqudeKNRLRiOffrrKozq30VLXt7Z/jhJdB92v7QTEog+4Uh+7oWVLKD0kpz0gpz8zPz+/qg155ZIKffdsrOgcSO+FWy63kZSmsmLM8//5qGRXF8FvYYtx35yEdVsrC30v0r+X27JQe8vvVQVs4HaO3VqwxORZ3rYM2ZFfplM6H7c6PJ2maK4Hd8FvUZb+vU8pFXXc8nIHO5HCvDKopI7+Ui5NRsQztIHTHqfxRZVJXxSCCkYpy67de6S6V3S7dSk/ns0kK1abnGpKt7QbRiPA8oMJpUldKqaQ7g2JQBv0qcNz2+zHg+oDuHRinWm4ppWLY7O8l5sp14lHRqYO2UFlg4TfTD86VIioeLqgddO2UR5waixOLCM+w2W+looXTHIBf/h+6e5R7tb1btcy8mR/2OvS3sN0kEY2QirurvVMtd63ZYrPSUPBw1eSfHIsT65vHmR9P+h5YrKI7jgY9gO74Tep6RndeuqMwmIKRitqpO1XmFAdTv34LDgZdYWGaVa7rNah0Duq2yb+13aDRkiNn0B8CTgshTgkhEsD3AJ8a0L13RX/YX9g29uLw65QqlRb50s6JIesz3XYa7MjhsWOb/T6wu06pcgxdzmHFZiQifA/XUPWynPblVgk7Y9EI4z4HXbvtGui106Bd/okx705p3dveKdfNe6qm67xSXm4nTs26LEqxU/DZ8gL25gz46U7DXPzUP5im4lHSiajnaUkqNfTgvG2Hiu6oHHSd7xxj1yu/35bb4F+hA8Y82uRYfFf9dlAMxKBLKZvA+4DPAU8DfyKlfGoQ994t/SGscspCcaR3qkeeGU9QbbTZ9lgYs7XtfMByv+zQ6zGsFWtEhPv2pxZWLbRbPXml3qTaaDsucjAOGXZXapWJIUv+RkuyaatWUQ07/Sbm3A7b7VQ2+cjvJzvAXL/uqBpEhb108iUX3ckk2Ki4n0lqrZD2M4huujOejJH2KJUF/zp66/t0Kn/0OyQ5iO5sbTd6dn1cV9Ad66Brz35bcnEGFHVHpY58h90p+lfoDJKB1aFLKT8jpbxdSnmrlPIXB3Xf3bI4keqsiITuDLyyl+UZNjvX9Ha9I2cv3euAZTtGaaVg2bb6bq1UY3Y86VvLOpHyPug65zKhC9673YG6l2VNWPfIr5ByASPt4hX2uy0OU/KyFOYvAA5NJDtL3CFYysL6HDfcVmzOZIxd/9wMqtcBy3as1csrfbqj0vZWyavboOK1uMpr61sIrjvWitRmq02uXFd2Bvzavr+G3pIdfKI7Zd1J7ei3MGIe+s3I4kSK1WK1o5yqI6XKJkUblZ1hJ3QXi7gZFZXlw2B4GwvZVE+nXC2ozZT7bYNqeVFOi6t8DbplVHwU21pabclfrhkra5U7pWfYXGcsHmWsb3JqaiyOEO4HDIBaDtqQv9cZCJKyAJ+w36HkEvxPS1JZ2AJG+mNyLN6nO1Vl3fE66DrnkrKwnvM63NzSHb+yy4U+3cmXzZW1KgOSz0HXGxX3VKn1WW6oeugLE72b3FmPtUHfI4tm2G8ZsG7plnc9aCYRJRrxXm2ZKznvGjhjLhZx65Qqq1wt+hVD2cvyWf7vtSeJUZ3iPbHldsCyncU+L6sTHSkYFb9Nltw83JhD/rIfv82hLBYnUj07Z1qTlSrpLnBPuRiHZDjLP+1jVLwOWN4pf3J3HrrP8n+nGnoLv/mXrW33A5Z7ZbcijOC645euy5Wc5y/GkzHiUeEdYXic5WpnhyNZMs5N6C+g2C8OrEE/NNkb9q8VayRiEd8vxVht6R72N1ptCtWmYw66s4eMj5e1q9BNMWXhZ1TcaujBMOiFatO1jtrrgGU7lpdlb3v7837y+xl0tyXsMwpzAH7pLtiZtlgrGpuiWQdSuNE56NpFdwrbTZoONfRgC/v9dEdxQFruizAC6Y7LSl1Lr53O5PWbf1EpuYSd6boguuO3StrNGRBCeJ5nCmrlxmDI32jJTluolFwOkgNr0K2VWdaxaqvmTLlKw3qN9G6TcuAfuqksH7ZYnEixvFVFSmlbuq2m1PbP6sdrx8nOKkmXjqmaskjGokyn451jySwPVzWP61f257ZB1qyHl2itclVLuRhy3ujTHT+sg67d5HeroQcV3VFL14GlO8ZRa9v1FsVqU7ntwT+6c0s3VhttKi7pGtWUxVQ6TiIW6cjf0Z0BeOheWwjP9FU22ak2WtSbbcV+azozNt0Pa0IUDrBBP26el3nFPLz2cr7SOUfQD68Vc9bI67Rz3UQqRiziHrp1F+b4e4nHpsco1ZpsVhosF6o0WpJj5t/kha+HXnGuoQe1sF/FSwHj8OCrG922B5Tl9zroOu+xa6DTsnGL7UaLZlsqdcp+3bkSSHfcyy478xeOBlEx5aKgO8en06wUalQbLa5sWG3vL79fhZdVQ++UNumkG10GVFXdEUJwbHqMK3nDoF/OV4hGRCfi9pPfd/7CRXdmMnHXOYAgqdJjjnbHX+8HxYE16HPjCdKJaOc08ku57knlfnh5iZYH6Ba6TXuEbkHC5pOzGUPufKV7SryC/H6rFd1q6MG/llvVywLjNHXrJPvLuQqzmQTjCnlElTyuq4c+7l5poTqpCMbh09GI4HK+gpSSy/kKJwLojlvbu5XNQbeW2y/loiK/pSdXNyqd78DSJy/8qnScdiq08CsI8DpguZ+TM+mefnt0asw39w7eB103W202K+7nCbvtUgnqFTrQbftL+QqNVpvrm1VluzMIDqxBF0JwYibN5VyFcq3JeqmmpNTgncf1O2zXq3xLZbc/i45i5MpcNjvliRl/xcj6HHTtddiuZWjccqF+ByzbOTmb5trmNo1Wm0u5YAbR+Kyd7b9db7HdaHmEzc6HU1iy2+/vRTwa4chUiku5CmulGpV6i5MKbQ/euuN3wIpXlVFh2/uAZTsnOrpjcwYU5Pero/caTDvOgEe6Tl13MlzOlTuDqapB9Krw2vCooQf3k6qg25dU5l+yqTgzmQSXchWubWzTakulfjsoDqxBB9NLzFc6Ib9qw3qVP/kddODXKb0OWLZjhf2XcxUu5SvEIoIjU/5hc9znoGu3GnpLduM97l6iasrl5EyGVltyfXPb6JQBDCI453H9Dtvt1HI7dOggHq4l/6V8pTOYqjoDXgek+B2w4uUMWNGRyhyQ1daXcobuZ1Mx18Oh7WRT3gddO20ZYZcd3KM71TpuMPppud4iV64bzoByv3XXHZV+61YQEMRDt+S/nC93ogxV3RkEB9qg376Y5cJ6mYcvbXR+V8F7YstQjCmXL3fay6ArLB+2GEtEOTGT5tyNAueuF7htYVx5g3wv+b1OPbI6vVceVFWpTy+OA/DACzmubW5zOkDbg7OX6HfqUbdsdGcu1O8s135OL47z7HKBx69udX5Xweuga7caeotpj1pu1ZJLMIzTbCbR0Z3bF7NKA0EkIoytFzx0320w9Zp/UTlg2Y7VT//m6RW2thuB+i04606nht5t/sVl22QIlkMHuH1xnHPXCzx13dSdBTXdGQQH2qDfc2KKVlvyh1+7xHgyxm2KDTuRilNvth13X9so15lK79xcycI7dFMrm7PL//ClDR69ssk9J6aUr/ML+906pdNeFBYqByzbeeWRCRLRCL//wEUAZfm9wma/w3a98rhBSkYB7jkxTbXR5o8fusJ8NslRhejIun9h23nrBbcadAujdM75ewuiO0II7jkxxYMX8jxxbYt7jk8pXQfuuuNVQw9GOsKtllvlgGU7dx2fRAj4va9eBNR1xytd1zkpyiPlAs5rSHajOxuVBv/1keucmst4npA2aA60Qb/bVORnlovcdWxS2cP1mhzyO2x3JpMwd1hzCt3UvRSAe45PsVqssbXd6PwtKrhN6nrV0Fu41RMHDTuTsSh3HJ3gmeUiAHcdm1K6zmvHQr/Ddr1quVU3FrOwjOCzK0XuPj6lXEc8ORan2ZZUHPbzyXnU0IOVcnHx0BXruC3uOTHN5XyFWrPN3QGcATfdKVTda+jBLAhwOdw8SMmlJcPphXGeWS6SiEV4+aEJpeu8KrzcNuaymPHy0ANGd9YAZOlOmBxogz47nuRbXr4AwA+86ZTydV65uFy55nnYrlctd5CUC8A77+6eEfK2Vx1Wvs4tj9vdD9p7k36nPKjKAcv9WG3+zlcfUapwAe8cet4nB+01B2AZFb9VrhbHZ9LctzQDwA+8cUnpGvDP43odku1Vyx1Ud9597zEA4lHBt7x8Ufk6t60XVA7JdtopEdQOWO7nn7/R0J1/fOaY74IuC68KKa8aevD20AvVJql4xHMvfzsvW8zy8kNGmuj733BS6ZpBEc561CHyG//zvTx1fYszZudUwWtP9Hy5zqk590mObi6useMMx63tBksBJkhmMgn+9qfeYmwopBjugaHYT99wGoyslIW7hz6dSXQmAu0EDTsBvuOuwxyaSHHXsUnla8biUWIuWy/kynViDpsrWcz4hM2ZRNQ1VebE7/3Aa3l2pci9J6aVr7FXihye7E3T5Mt1z1y8NdAaZ772/o1BJhXBWCn9xZ94M5lkzDVn78TEWIyL6zu//3xnUZS77rjt5xKkXNfie+87zq3zGe7ZRdu7RdZuNfTgPQewVQnW9kII/uS9b+DCWplXaw99sIwlooGMOXjncXOlumfKomtUdiq2Xw7ViaW5DMcDlj255UFVvCy3SguvTb3cEEJw36kZpaoe+zVu8udKxhJ8t/SHVy33ZsW9QsONTDIWyJiDu+5Iaaz2dZu/APc5ACklG5WGZ2TlxG0L2R2Dih9uqy3XPbaMsHCr8Oqc8BVQd153y6yydw7eB1277eNiYazNcHYGnM6w9WMiFQ/dmMNLwKDvBrfZcuuA4jnPPKhzp2ya+WuV8rG94nbQ9bq58ZaX/G613Bs+IesgcavSMfZx8V5G7bYnR77ine4YFG5zAJV6i1qz7evhwk6jUjDPUQ1FfpeDrv0mpMG9IMD6PkLRfRfdyXmU6wJEI4Kpsbij7mxU6qHIPgi0QXfALY+7ud2gLf3ziLDToFtnXQYd6XeD20HXfpOKxmsJx1pua2GG1/zBoHCbA8iVvQdT8IowGqFUG7hNqKtGR7CzlnvTY8uAQeN20LVaDj3pWMtt6U5Y8rsNSH59zz3CcF9herOhDboDbrk4a2tZLy/LGsn7FSNsDxd2hv25Up2IcK+hB/eU0YZ5OIDqpOJeMHa7dA6b/TqW2zqAjXKdmZCiI9jpDKhER535l74JdRVjOijcto5YL9XIJmOeE4NWSqjfy82X62RTMaXl+3vFbafUXEktunNLlYbRbweBNugOJGIRxuLRHaFb9zxL9y/XrZY7zE7pVr5lLfv32v7WLcIwUhZx361zB4HXHIDfnuSuXpbHKsdB4nbQtUp05FbLvZv5i93ipjteOxVaWH+bk/xhebhOcwAtM1Xq1W/BWXesNKv20EecibGdI31nc6VdhP1Wpwwnh+58yEW+XPM1iG5zAMaCqnCU2mnXvGqjRanWVGr7ftnrzTbFWjOUdFHnoGs33fEwDFYtd3/KJV+2UhZh5KCdS//8JhXBvZY7P2Td2azUaUvvtgdrg67eawtVI80alvx7RRt0F5wWWPgtTrBwmpizFCXMsHmHh66QsugcmOvkZYXVKc22t0/M5hWiIzA65Xaj1XNQ92aIHi4YA6prdOfr5e5c2NU9WCKcSVFwi+58nIFx50ndzUojlHQXuPVba1M0P2fGKAiwn6najaz0pOhI4zS5YpVu+eXTnPK4XnthDxqviTnfsDnt3Ck3yg3Hk2r2g8kx46DraqM7uZbz2LbYjtN+LvldlM3tBScvMVeqkYpHdtSX9+MU9ucrRv19GMeYuemO27GLdizddko3hjWYWuk6uzOgUnIJRr/tP6g7zH47CLRBd2EqHe/MzlvkzX1c/CZ3HFMu5TrpRDRQTfZusdI6/fKvl2rM+Sj1WCLKWHxnLbfb4cb7QVf+rgw5hYUt4FzLbe3jEVbp2VQ63okKLFTy/+Bs0K0a+jCOMZs028guv5TS8+g/Cysl5JhDD8kgTqWNg67tKSOVkktwXi26EWJkPQi0QXdhNpPcseotV64pfbFWyqUnZRBSHTRAOhEjFY/0yF9v+u/jYtGfMrI2ZgpLfqc9WVRy0OA8qbubhS17YXZ853Fm6woGEUxnoO+gbq9TmgaNU4TmdRaqnVg0YjhCtmurjRaVeis0D33WIWWYC5Aqhd45gLz20A8GM+MJcqVeo5wr1ZlTNIjNtuzxEowqi/DycLOZZM82uJZRUzIqfSf/lGpGhw7NoJsyrtsMW1Avy96h/faAGTSGUe730P1TFuBcy71RboQWXVhG2S6/ZRDnFM717I8wNkOsQTc+36y0selOJ13nI4PTwq4w5y8GgTboLsxaRtlWreB12o8dJy8xX2mEOsrPjidY34WHCzvP5rTCztC8rE6n7MqwXq6RMCtIvHDak6O7UjE8+ft33PTbMsKiU8td6fUSwwz5+3d97BzMoaL76d5rw55UtPR7vW9A8try2sK539ZJmIfGjALaoLtgeSP9yqni4XYVo3vtbvZx2QuzfRsldT1cf6PSX/rXnVQMqVM6hM35ktH2fnlkp1rufKVONhkLtC/IXrDktwYSKaXSKldwmwMIb1IRzJRRycEZUNR9p3RXWM6M1W/7IzS16MjZGZjOqJ0UdTMwEA0XQvyKEOIZIcTjQog/F0JMDeK+w6Q//Ko32+TLdaWw03rPWrHXoIfrofd2ytVi1ZTNX4a5bJK1Uq2TbgpzlSvAeDJGIhph3TYgrRZrSm0vhGA2k9zZ9iEPptD1EgvbTerNtqLuGNda8ncWtoSpO32T+mum7syryJ/tbfswF9QBnbSmPeWyWlDTnVQ8SjYZ65M/3Mh6rwzKZfkCcKeU8i7gOeBnB3TfodGpqTUVwzKIhyZTrtdYWO9Z3jKuabSMhS2hGvRM7xzAckFd/sWJFPVmu1MlE7aXJYQw8vi2AWmlUGVxwl92gMXJFCvm3wvmPi4hbq402xfdWW2/uAvdKVab5sKWMOXvnZhdLlSJRoRSdHd4IsVGpdE57WuzEm66KxmLkk3Fegak5UKVwwptD8Z3ZLU9mBVGLzWDLqX8vJTSSjZ/DTg2iPsOk/k+L9syEIcUjMpMOkE8KlguGNcGCVkHxXw2Sb3V7swBrGxVyaZivnXQ0P0bLcW22iBM+efGez295UKVQ5P+BgXg0ESyp1Oul2pKxmhQ9HvZywF0xxq0LH1bK6lPSA6KufEkG5UG9aYxB7C8VWMhm1Q68csatFZN3V8r1oiIcFa5WszbdEdKyWqhpjSYgvEdLRf6declZtD7+BfAZ51eEEK8RwhxVghxdm1tbR8+enDMjieJRQQ3TMOwvGUoiIqXGIkIFrJdL/HG1jYAR6bUlGoQdDw9U4blQlXJoBjXGsajK3+V8WSMbIBN/vfKIZunVG202Kw01OXv65TLW1WlyGRQWHuQW7qzsqVu0FPxKFPpePd7M69V9TAHgfVZVlQaJDrqOAM23ZnPJgMdLLJXDk2mOn0uX65Tb7WVdWdxottvpZTc2KpyRPE82ZsB5VYWQnxRCPGkw7932d7zc0AT+KjTPaSUH5JSnpFSnpmfn9+79PtINCJYnOgalSApC+t9nWs7HTo8xbA65XVTsZcLNWXZF/s65fKWesg6KI7YOqXVwYKkXIrVJpV6k2qjRa5c50iI8o8lTKO81TWIAAsTqhFGquNAWG0Q9KCKvdA/IAVzBnY6EmHKDob8O/ptAGdmtVij1ZZsVhrUmuqDwc2A8lpiKeX/5PW6EOKfAe8AvkU6HXk+ghyeTNk89G0SsYhy6HhoMsVT17YAuD4EL+uQ2YmWbfKfXlAbRBeyKYSAG5uGMblRCNfDBUP+QrVJudbsfAeqMnQGs80q8ajo3C9MDk10dedGocp0Oq68StjuYVr3UB0MBoHVfje2qkgpWd6q8k23zSlda31HHd3ZqnJ6wf3Yvf3g8GSKFdMoW/qvnHKZHKPVlqwVa505kLCdmb0wqCqX+4GfAd4ppawM4p43A/aOdSlX4cRMWrl86cRMmmub2zRbbZa3tknGIqFObC1kk4ZR3qqyXW+xUqhxUvEou0QswpHJMS7nja9yeWs7dKW2GxXrjNOTM2rnsZ4w/84r+UrHIA5Dfkt3LucqnAhwluyJmTSX85VOyD+bSYSyZYRFd2J2m3y5TqnW7LSpHxOpOFPpuE13huEMpDpG+VJHd9Tkt/7Oy/lKN7J+qRl04DeALPAFIcSjQogPDui+Q+XwZIrrW1XabcmlXIWlWTWlADg1m6HRklzfrHJ908jDhVnLGo9GWMgmub65zaV8GYCTHodb97M0l+ZirkKt2WK1WAvfw+142dtczJWJRYTyHIR1EPeF9TLXTU9xGBGG9dkX1suBdGdpNkOx2iRfrnN9c5vDIc69AGRTccaTMa5vVrmYM3RnaS6Y/JdyFbYqDUq15tCcgWub21zKlckmY8plk6dM3blo051RyqEPZPs2KeVtg7jPzcYt8+PUm22ubFS4mCvz906rhZ0AJ80OfCFX5vxqiVsCGNNBccvcOOdXS1xcNzrlqQBe4snZDJ954gYX1stICbfOhyv/LebnnV8tcTFX5sRMWnlibSaTIJuMcTFXZrwUIxYRHJ9WN0iD4Nb5DBuVBje2trm+tc3SrHrhl2U8L5q685qTwQ6qHgS3zGc4v1riwrrh4S4F0J2l2TQPXdzg/FrRuNdcuCmXW+aNz3thtcSFXIWluYyyM3VkKkUsIriQK1OpNRlPxljIhpfu2it6pagHty9mAfjyc2vUmu1AHu4p873PrxR5cb3EafNeYXL7omHQX1izPPRgEcZmpcGDF/LmvcKVf348yXQ6zvOrRV5cK3cGSBWEECzNZbiwXua5lRKn5jKhrRK1sNrri+dWkDK4hwvw5LUC1za3uX0xXIMIcHohy3MrRS6sl4gIOBZgQFyay3B9a5snrhpzSC87FK7unJhJk4xFOvIH0Z1YNMKJmTQX1gzduW1hfGRWiYI26J6cNjvSxx68AsBdRyeVr53PJlnIJvnkN67RaMnhdMrFLKVak796/AZLs+nO4QUq3Gn+rX/80BWiEdHxmMNCCMHpxSzfuLTJcytFXhWg7QHuPDrBo5c3eWa5EPpgBF2DbunOq45OKV97cjZDNhXjjx8yrh2WM7BarPGlZ9Z4xeGJQAPinUcmkRL+9OGrjMWjHA05ZRGNCG5bGOeBF3JcyW/vQncm+cblDZ5fLQ6l3+4FbdA9mEjFuXU+w7kbBcaTMe44MqF8rRCC190yy7kbBQBefXxqn6R0527zM8/dKPC6U7OBrr3nxBSJaISnrhd4xeGs5+HA+8U9x6d4dqVIW8Lrbgkm/+tOzVKsNbm6sc2rjwfr0INgcSLJoYkU524UmBtPBkpZRSOC1y7NdHTnrmPhy78X3XntqRmEgKeuF7jr2GQo59D2c/fxqU77BdadW2ZYLdZYL9WH0m/3gjboPnzXvUbu8++/fCHw4oi33XkIgNML49w6H/5If8eRiY539NY7FwNdm4pH+SZzzuC77hnOwl+r7QHuPREsj/xNtvmO77z76MBkUkUIwXfeY3zuW+9YDBy232/qzmtOTodexw3w2qWZzoHX33ZHMN2ZHItzj2kI333vcHUnEY1wZwBHDOAtL1voPH7HXUcGKtd+I4ZVMn7mzBl59uzZoXx2EKqNFl84t8K33bG4Ky/1gRfWOTI5xtIQJkUBzq8WyZcb3HdqJvC11t/+1jsOhZ6DtvjSs6u8bDG7q0qDi+tlrm9t88Zb1SezB0m51uRLz65y/x2HdrVS8svPrXHLfCZQ/nqQPH2jwHajFXgwhb3/7XtFSsnfPL3KXccnWcgGr7LZS7/Zb4QQD0spzzi+pg26RqPRjA5eBl2nXDQajeaAoA26RqPRHBCGlnIRQqwBl3Z5+RywPkBxwkbLPzxGWXbQ8g+Tm0X2k1JKx42ZhmbQ94IQ4qxbDmkU0PIPj1GWHbT8w2QUZNcpF41GozkgaIOu0Wg0B4RRNegfGrYAe0TLPzxGWXbQ8g+Tm172kcyhazQajWYno+qhazQajaYPbdA1Go3mgHDTGXQhxP1CiGeFEOeFEO93eF0IIf6j+frjQoh7Va8Ngz3K/7tCiFUhxJPhSt35/F3JLoQ4LoT4khDiaSHEU0KIHw1f+j3JnxJCPCiEeMyU/xdGRXbb61EhxCNCiE+HJ3XP5+9F7y8KIZ4wTzsbyn4ge5R/SgjxZ0KIZ8w+8IZwpbchpbxp/gFR4AXgFiABPAa8su89bwc+Cwjg9cDXVa+9meU3X3szcC/w5Ii1/WHgXvNxFnhulNre/H3cfBwHvg68fhRkt73+E8AfAZ8eJd0xX7sIzIUt9wDl/33gB83HCWBqWH/Lzeah3wecl1K+KKWsAx8H3tX3nncBH5EGXwOmhBCHFa/db/YiP1LKLwP5UCXusmvZpZQ3pJTfAJBSFoGngbD3rN2L/FJKWTLfEzf/hVktsCe9EUIcA74d+HCIMtvZk/w3AbuWXwgxgeGI/WcAKWVdSrkZouw93GwG/Shwxfb7VXYaBrf3qFy73+xF/mEzENmFEEvAPRhebpjsSX4zZfEosAp8QUoZpvx7bftfA34aaO+TfH7sVX4JfF4I8bAQ4j37JqU7e5H/FmAN+D0z5fVhIcRw9srm5jPoTqcA9HtKbu9RuXa/2Yv8w2bPsgshxoFPAD8mpSwMUDYV9iS/lLIlpbwbOAbcJ4S4c7DiebJr2YUQ7wBWpZQPD14sZfaqO2+SUt4LvA34YSHEmwcpnAJ7kT+GkSb9bSnlPUAZGMr8Hdx8Bv0qcNz2+zHguuJ7VK7db/Yi/7DZk+xCiDiGMf+olPKT+yinGwNpezNc/lvg/oFL6M5eZH8T8E4hxEWMVME/EEL84f6J6sie2l5Kaf1cBf4cIwUSJnu1O1dtEd2fYRj44TCs5L3TP4zR7kXgFN3JiTv63vPt9E5OPKh67c0sv+31JYYzKbqXthfAR4BfG1HdmcecyALGgK8A7xgF2fve8xaGMym6l7bPAFnb4weA+0dFfvO1rwAvMx9/APiVsL+DjizD+mCPxn07RpXEC8DPmc+9F3iv+VgAv2m+/gRwxuvaEZP/Y8ANoIEx8v/LUZAd+CaM8PNx4FHz39tHpe2Bu4BHTPmfBH5+VGTvu8dbGIJB32Pb34JhQB8DnhrRfns3cNbUn78ApofxN0gp9dJ/jUajOSjcbDl0jUaj0ewSbdA1Go3mgKANukaj0RwQYsP64Lm5Obm0tDSsj9doNJqR5OGHH16XLmeKDs2gLy0tcfbsUPbh0Wg0mpFFCHHJ7TWdctFoNJoDgjboNtZLNdZLtWGLodFoNLtCG3QbP/vJJ/jpP3t82GJoNBrNrhhaDv1mZLNSp9nWC600Gs1ooj10G822pKUNukajGVG0QbfRaksaLW3QNRrNaKINuo1mS9JqD+uMAI1Go9kb2qDbaLUlTe2hazSaEUUbdBvNdpuG9tA1Gs2Iog26jVZb0tIeukajGVG0QbfRbEsauspFo9GMKNqg22jpskWNRjPCaINuo9mWNFo6h67RaEYTbdBtaA9do9GMMtqg22i22rpsUaPRjCzaoNtotaUuW9RoNCOLNug2mm2JlNDWaReNRjOCaINuw8qfay9do9GMItqgm0gpO1vn6olRjUYzimiDbmK34XrHRY1GM4pog27StKVZmroWXaPRjCDaoJvY0yw65aLRaEYRbdBN7EfP6f1cNBrNKKINuol9l0W946JGoxlFtEE36fXQdQ5do9GMHtqgm9jz5nr5v0ajGUW0QTfpqXLRHrpGoxlBtEE30R66RqMZdbRBN7Hn0Ju6ykWj0Ywg2qCb9HroOuWi0WhGD23QTexpFr2wSKPRjCLaoJu09MIijUYz4miDbqL3ctFoNKOONugmLT0pqtFoRhxt0E2aumxRo9GMONqgm/R66DrlotFoRg9t0E20h67RaEYdbdBNWnrpv0ajGXG0QTexe+V6UlSj0Ywi2qCb6L1cNBrNqKMNuoney0Wj0Yw6SgZdCHG/EOJZIcR5IcT7HV5/ixBiSwjxqPnv5wcv6v6i93LRaDSjTszvDUKIKPCbwLcCV4GHhBCfklKe63vrV6SU79gHGUNBe+gajWbUUfHQ7wPOSylflFLWgY8D79pfscKnp8pF59A1Gs0IomLQjwJXbL9fNZ/r5w1CiMeEEJ8VQtzhdCMhxHuEEGeFEGfX1tZ2Ie7+0dQLizQazYijYtCFw3P9Luw3gJNSylcD/x/wF043klJ+SEp5Rkp5Zn5+PpCg+43ey0Wj0Yw6Kgb9KnDc9vsx4Lr9DVLKgpSyZD7+DBAXQswNTMoQ6KlD15OiGo1mBFEx6A8Bp4UQp4QQCeB7gE/Z3yCEOCSEEObj+8z75gYt7H5ieehCaA9do9GMJr5VLlLKphDifcDngCjwu1LKp4QQ7zVf/yDw3cAPCSGawDbwPVLKkbKKlhFPxaJ6UlSj0YwkvgYdOmmUz/Q990Hb498AfmOwooWLVeWSikf0pKhGoxlJ9EpRE8tDT2oPXaPRjCjaoJu02pKIgHhM6By6RqMZSQ6cQa82Wru6rtmWxCIR4pGINugajWYkOVAG/epGhVd94HM8emUz8LWttiQaEUQjQpctajSakeRAGfRrG9s0WpLzq6XA1zZbklhEEItqD12j0YwmB8qgV+pGumWzUg98bavdJhoVxLSHrtFoRpQDadA3dmHQjRy6IBbVk6IajWY0OVAGvVxvArBRaQS+1sqhGx66NugajWb0OFAGvVIzDPpuUi5WlUssohcWaTSa0eRgGXSzZHGjvAcPXadcNBrNiHKwDHptADl0nXLRaDQjysEy6HuYFG2126aHHqGhq1w0Gs0IcsAMendSNOhmj81Wd1K0pVMuGo1mBDlgBt3w0OvNNtsBtwBotSWxqF5YpNFoRpcDZtCbncdBSxebbUk0EjFy6LrKRaPRjCAHzKB3vfKNcrA8ektPimo0mhHnQBn0cr1FKm78SZuBPfTupKhOuWg0mlHkQBn0Sq3J0akxIHilS6+HrlMuGo1m9DhYBr3e4uh0Ggi+WrRpW1jU0CkXjUYzghwwg97kyGQKgHzA1aKWhz41lqBUa+padI1GM3IcMIPeYnIszngyFjjlYtShR5jPJgHIlYIvTtJoNJphcmAMeqstqTXbpBMxptLxwCkXy0O3DPpasbYfYmo0Gs2+cWAMulWDnk5EmU4ndlGHbhxw0THoperAZdRoNJr95AAZdKMGPZ2M7tpDjwrtoWs0mtHl4Bn0XXvoRsplbjwBaIOu0WhGjwNj0Ms1K+USYyaT2FUdejQiSMaiTI7FtUHXaDQjx4Ex6NZmXOmEkXIpVpuBFgg1zc25AOazSdZK2qBrNJrR4sAYdLuHPp020iab2+ppF8tDB5gfT2oPXaPRjBwHxqDbc+hT6TigvkGXlJLteotENAqYHro26BqNZsQ4MAa966FHOx666sToja0q240Wp+YzgDboGo1mNDkwBv3Ja1uMxaMcmkzZDLqah/7sShGA2xfGAcOgl+utziCh0Wg0o8CBMehffn6d198yQzLWTbmo1qI/bxn0xSwAhyaM/WAu5sr7IKlGo9HsDwfCoF/JV7iwXubvnZ4HYDoTLOXy7HKJ+Wyyc90bb5tFCPjCuZX9EVij0Wj2gQNh0L/8/BoAb759DoBMIko8KpRTLs+vFnmZ6Z0DLGRTnDk5zV8/uTx4YTUajWafGHmDXqw2+OB/f4Fb5jPcOm/kwIUQTKcTbCpsodtuS55fKXF6cbzn+fvvPMwzy0VeXCvti9wajUYzaEbaoBerDX7iTx7j2sY2v/zuuxBCdF5bmEhyKe+fA3/y+hbbjRavODTR8/y3v+owyViEX/yrp5FSH3ih0WhufmLDFiAoL66V+OLTK1zMVfj8U8vky3X+zdtfwZmlmZ73ffPt83zwv79IvlxnxsyNO/FfvnqRdCLKW+881PP8ockU//qtL+Pf/tXTfPgrF/hf33zLvvw9Go1GMyhGzkN/bqXEv/vMM3zq0eucOTnDn773jfzg39tpbN9252FabckXzrnnwVcKVf7y8ev84zPHmRyL73j9X7zpFG+9Y5Ff/MzT/OwnH1deqKTRaDTDYOQ89Le8bJ4nPvBtjCdjPSmWfu44MsGJmTR/+dgN/slrT+x4vdlq8+N//ChCCH7gTUuO94hEBL/1fa/hl/76GT78lRf51KPXeefdR/iWly/yhltnySRHrvk0Gs0BRskiCSHuB34diAIfllL+332vC/P1twMV4J9LKb8xYFkBSMWjpOJR3/cJIfgnrz3Or3zuWT79+HXecdeRzmvVRouf/rPHeeCFHP/hH72ak7MZ1/tEI4J/8/ZX8O57j/E7X3mRv3jkOh978ArxqOC2hSynF8a5fXGcU3PjLE4kWcimmM8mGUv4y6jRaDSDRPhN+AkhosBzwLcCV4GHgO+VUp6zveftwP+OYdBfB/y6lPJ1Xvc9c+aMPHv27N6k96HRavOPPvh3vLBa4j99/2t4461zPHF1i5/5xOOcu1HgZ+5/OT/0llsD3bPWbPHwxQ2+cn6dZ24UeG6lxLXN7R3vyyZjTGcSjCdjjKdiZM2f1u+ZRIxELEIyFiEZi9oeR8zHUZLxCImo8Xs0IohFBBEhiEWF+Xv3eftPr8hFo9GMNkKIh6WUZ5xeU/HQ7wPOSylfNG/2ceBdwDnbe94FfEQao8PXhBBTQojDUsobe5R9T8SjEX7z++7ln/3ug/zTD3+dpdkML66XmRtP8jvff4ZvfeVi4HsmY1HeeNscb7xtrvNcudbkUq7CWqnGaqHKarHGWrHGZqVOqdakWG2yXKhSWmtSqjYp1prUm+pb+wYlIug19lFBVAgiEUFEgMD8KQSRSPf3iBCIHT+t9/b+HhECQe97rXt1frc+QwAYg4wQxiNrzLHe331NWG813ydsj3uf7wxb1nUK98f+u9v7bJ/bI5fT/T3k6MdtmHV+r/O73cZqx6dd3rxfcgS5r/G8uuOxX23k/t5g9w7KXcemuO/UjP8bA6Ji0I8CV2y/X8Xwwv3ecxQYqkEHODo1xp//qzfy4a9c4LGrm7z7Ncf4p68/6TgJulsyyRivPDLh/0YbzVabeqtNrdH9WWu2qDXb5r8W9WabaqNNqy1ptq2fkpbtn/F72/jZMn5vy+77mq3u621p7CzZlhIp6fm9LUGC+Zqk3QaJ/Rps1+382WpLGi3z3vTd2xy7JPSUgEppfEb3cfd12flv5/M7rrMFmVJar7jdv/POzmPH+3vI4XV/bO+3I52epFd2j8s7f5vj84r31dw8vPebbx2aQXcak/rVReU9CCHeA7wH4MSJnROV+0U2FefHv/X20D5PhVg0QiwaIe1eUanRhIbTYOE2KDgPIOqDjdu9gwx6brjLvLe/z3i/wz0U5eonEd2fAkMVg34VOG77/RhwfRfvQUr5IeBDYOTQA0mq0Wj2Daf0R7D0gp63uRlQGSYeAk4LIU4JIRLA9wCf6nvPp4DvFwavB7aGnT/XaDSalxq+HrqUsimEeB/wOYyyxd+VUj4lhHiv+foHgc9gVLicxyhb/IH9E1mj0Wg0TviWLe7bBwuxBlza5eVzwPoAxRkUN6tccPPKpuUKhpYrGAdRrpNSynmnF4Zm0PeCEOKsWx3mMLlZ5YKbVzYtVzC0XMF4qck1cnu5aDQajcYZbdA1Go3mgDCqBv1DwxbAhZtVLrh5ZdNyBUPLFYyXlFwjmUPXaDQazU5G1UPXaDQaTR/aoGs0Gs0BYeQMuhDifiHEs0KI80KI9w/h8y8KIZ4QQjwqhDhrPjcjhPiCEOJ58+e07f0/a8r6rBDirQOU43eFEKtCiCdtzwWWQwjxGvPvOS+E+I9ij3vvusj1ASHENbPNHjW3Ww5bruNCiC8JIZ4WQjwlhPhR8/mhtpmHXENtMyFESgjxoBDiMVOuXzCfH3Z7ucl1M+hYVAjxiBDi0+bv4beVNHfXG4V/GCtVXwBuARLAY8ArQ5bhIjDX99wvA+83H78f+CXz8StNGZPAKVP26IDkeDNwL/DkXuQAHgTegLEZx2eBt+2DXB8AfsrhvWHKdRi413ycxdjj/5XDbjMPuYbaZuY9xs3HceDrwOtvgvZyk+tm0LGfAP4I+PSw+uOoeeidvdmllHXA2pt92LwL+H3z8e8D32l7/uNSypqU8gLG1gj3DeIDpZRfBvJ7kUMIcRiYkFL+nTS06SO2awYplxthynVDmqdoSSmLwNMYWzwPtc085HIjLLmklLJk/ho3/0mG315ucrkRilxCiGPAtwMf7vvsUNtq1Ay6277rYSKBzwshHhbGdsAAi9LcjMz8uWA+H7a8QeU4aj4OQ773CSEeF0ZKxgo9hyKXEGIJuAfDu7tp2qxPLhhym5kphEeBVeALUsqbor1c5ILhttevAT8N2E+uCb2tRs2gK+27vs+8SUp5L/A24IeFEG/2eO/NIC+4yxGWfL8N3ArcjXHoyf8zLLmEEOPAJ4Afk1IWvN4apmwOcg29zaSULSnl3RjbYd8nhLjT4+3Dlmto7SWEeAewKqV8WPWS/ZJp1Ay60r7r+4mU8rr5cxX4c4wUyooZLmH+XDXfHra8QeW4aj7eV/mklCtmJ2wDv0M37RSqXEKIOIbR/KiU8pPm00NvMye5bpY2M2XZBP4WuJ+boL2c5Bpye70JeKcQ4iJGGvgfCCH+kGG01V4mAcL+h7Hd74sYEwnWpOgdIX5+BsjaHj+AoeS/Qu/kxy+bj++gd/LjRQY0KWref4neycfAcmDsd/96upMwb98HuQ7bHv84Rv4wVLnM+3wE+LW+54faZh5yDbXNgHlgynw8BnwFeMdN0F5ucg1dx8x7voXupGjobTUQwxLmP4x915/DmBn+uZA/+xbzi3gMeMr6fGAW+BvgefPnjO2anzNlfZY9zqL3yfIxjNCygTGy/8vdyAGcAZ40X/sNzNXDA5brD4AngMcxDkM5PAS5vgkjfH0ceNT89/Zht5mHXENtM+Au4BHz858Efn63uh6SXEPXMfOeb6Fr0ENvK730X6PRaA4Io5ZD12g0Go0L2qBrNBrNAUEbdI1GozkgaIOu0Wg0BwRt0DUajeaAoA265iWDEOJHzF0NPyqE+E4hxM+bz39ACPFTfe+9KISY87jXF+2752k0NwOxYQug0YTIv8Ko+b0ghHgAeOce7vUH5v1+cSCSaTQDQHvompcEQogPYiwM+5QQ4meAmpRyXeG699r22L4ghPiS+dKngO/dR5E1msBog655SSClfC/Gvhh/H2NPjW/0veXHbYb7UeCIed0HpbER1GsxVr7+qvn8BpAUQsyG8xdoNP7olIvmpchhYK3vuf9XSvkfrF/MjZbs/Drw36SUf2l7bhXD8Of2Q0iNJijaoGteimwDk6pvFkL8c+Ak8L6+l1LmvTSamwKdctG8FHkauE3ljUKI1wA/BfxTaWzNaj0vgEMYRxJqNDcF2qBrXop8GbhH8QDe9wEzwJfM/Lp1xNhrgK9JKZv7JaRGExS926LmJYkQ4teBv5RSfnEP139KSvk3g5VMo9k92kPXvFT5d0B6D9c/qY255mZDe+gajUZzQNAeukaj0RwQtEHXaDSaA4I26BqNRnNA0AZdo9FoDgjaoGs0Gs0B4f8HbS17/cO3qWYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#导入使用的库\n",
    "import numpy as np\n",
    "import pylab as pl\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "f = 8000\n",
    "N = 512#自己改N，相应长度也要改\n",
    "t = np.arange(0, 1, 1/f)\n",
    "x = np.sin(2*np.pi*100*t)  + 2*np.sin(2*np.pi*300*t)\n",
    "A = x[:N]\n",
    "y = np.fft.rfft(A,N)/(N)\n",
    "freqs = np.linspace(0, int(f/2), int(N/2+1))\n",
    "B = np.abs(y)\n",
    "plt.figure()\n",
    "plt.subplot(211)\n",
    "plt.plot(t[:N], A)\n",
    "plt.xlabel(\"t/s\")\n",
    "plt.subplot(212)\n",
    "plt.plot(freqs, B)\n",
    "plt.xlabel(\"f(Hz)\")\n",
    "plt.show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b65289e1",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
